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We explain why compact (7(1) confines and how to fix it. We show that plaquettes of negative trace carry most 
of the confinement signal in compact 517(2). We show how to perform noncompact gauge-invariant simulations 
without auxiliary fields. We suggest a way to simulate fermions without doublers. 



1. INTRODUCTION 

We reort here the results of four different stud- 
ies all of which are attempts to find better lattice 
actions. The first study explains why compact 
U(l) gauge theory displays confinement at strong 
coupling and shows how to remove this artifact. 
The second study shows that most of the confine- 
ment signal in compact SU(2) gauge theory is due 
to plaquettes of negative trace. These results ob- 
tain for both the Wilson action and the Manton 
action. The third study shows how to perform 
noncompact lattice simulations that are exactly 
gauge invariant and that do not involve auxiliary 
fields. The fourth study presents a lattice action 
for fermions that avoids doublers. 

2. WHY COMPACT 17(1) CONFINES 

Compact lattice simulations of f7(l) gauge the- 
ory display confinement at strong coupling, as 
shown by Figure 1 which plots Creutz ratios |Q 
obtained with the Wilson action. This lattice ar- 
tifact is obvious in the figure at (3 = 0.25 and at 
[3 = 0.5, and incipient at = 1. It arises because 
U(l) is a circle; if one cuts the circle, then there 
is no confinement, as shown in the figure by the 
X{i>j)' s labeled "cut O" which follow the curves 
of the exact Creutz ratios down to j3 = 0.25. 
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We ran on a 12 4 lattice and began all runs from 
a cold start in which all links were unity. We used 
a Metropolis algorithm and rejected any plaque- 
tte whose phase 9 was either greater than tt — e 
or less than — tt + e. We saw no confinement sig- 
nal as long as the step size was smaller than the 
thickness 2e of the cut in the phase 9. Simulations 
with Manton's action show similar results. 

One may interpret these results in terms of 
monopoles. Since the phase 9 of each plaquette 
is required to lie between 7r — e and — tt + e, it 
follows that for e > 0, no string can ever pen- 
etrate any plaquette. For the Wilson action we 
took .02 < e < 0.1 and noticed no sensitivity to 
e within that range. For the Manton action, we 
took e = 0.026 in all runs with cut circles. 

3. COMPACT SU{2) 

In view of these results for U(l), one might 
wonder whether similar lattice artifacts exist in 
the case of the group SU(2). Inasmuch as U(l) 
and SU(2) have different first homotopy groups 
(vri(J7(l)) = Z but 7ri(S*C7(2)) = 0), one might 
assume that excising a small cap around the an- 
tipodc (g = — 1) on the SU(2) group manifold 
(the three sphere S3 in four dimensions) would 
have little effect on Creutz ratios. 

To check this assumption, we ran from cold 
starts on an 8 4 lattice and used a Metropolis al- 
gorithm in which we rejected plaquettes that lay 
within a small cap around the antipode. The step 
size was small compared to the size of the ex- 
cluded cap. As shown in Figure 2, the Creutz ra- 
tios \(h j) do not depend upon whether the small 
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Figure 2. The SU (2) Creutz ratios \{h j) as given 
by Wilson's action with and without the exclusion 
of a small part of the SU (2) sphere. The curves 
represent the perturbative Creutz ratios, x(2, 2) 
(solid), x(3,3) (dashes), and x(4,4) (dots). The 
X(i, j)'s substantially overlap. 



Figure 1. The U{1) Creutz ratios x(hj) as given 
by Wilson's action and by Wilson's action on a 
cut circle. The curves represent the exact Creutz 
ratios, x(2, 2) (solid), x(3, 3) (dashes), and x(4, 4) 
(dots). At /3 = 1, the symbols for the Wilson 
action are plotted to the right of those for the cut 
circle; at (3 = 1.5 they overlap. Wilson's action 
on the full circle confines for (3 < 1. 



cap was excluded. 

But what about the excision of a large cap? 
To study this question, we again began with cold 
starts on an 8 4 lattice and employed a Metropolis 
algorithm with a small step size. We rejected all 
plaquettes that had a negative trace, thus exclud- 
ing half of the SU(2) sphere. 

In Figure 3 we plot the Creutz ratios x{hj) 
both for the usual Wilson action and for the 
positivc-plaquctte Wilson action. As shown in 
the figure, the x(*> j)' s of the positive-plaquctte 
simulations exhibit behavior that is substantially 



3 



10.0 F 



1.0 



0.01 



0.001 



; 1 




T 

"a © 


- T- A q - 


$ ——J 

r T i 


I A 

1 : 


"** _ < 








°x(2,2) 




= a X (2,3) 


5t/(2) 1 


I vx(2,4) 


- <>X(3,3) 


Creutz ratios 


" > X (3,4) 




■ <X(4,4) 
i 


i 



1.0 



1.5 





2.0 



Figure 3. The SU(2) Creutz ratios x(hj) as 
given by Wilson's action with and without the 
exclusion of half the SU (2) sphere. The curves 
represent the perturbative Creutz ratios, x(2, 2) 
(solid), x(3, 3) (dashes), and x(4,4) (dots). In 
the two cases, the x(hj)' s substantially differ. 



perturbative. The confinement signal has disap- 
peared. Apparantly the plaquettes of negative 
trace carry most of the confinement signal. 

4. NONCOMPACT SIMULATIONS 

Suppose we write the fermion field 
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(1) 



in terms of some orthonormal basis vectors e a {x) 
that may vary with position and time 



ip(x) = ip a (x)e a {x). 

Then the derivative d ll ij){x) has two terms: 
d^ipix) = e a (x)d f _ l ip a (x) + ip a (x)d fl e a (x). 
And so if we let the gauge field A a Mx) be 



(2) 
(3) 
(4) 



Af{x) = (-i/g)ei(x) ■ d^e b (x), 

then the free action density iip^^d^ip becomes 

C v = i&y {8 ab d^ + iAf) rp b = ii> a ^Df^ b (5) 

as in a gauge theory. 

Under a gauge transformation 

^a( x ) = 9ab(x)tpb(x), (6) 

the field ip( x ) ail d therefore the action is invariant 
if the vectors e a {x) transform as 



e'a(x) = 9ca (x)e c {x) 



(7) 



which for unitary groups is e' a (x) — g* c (x)e c (x). 

To generate U(l) gauge fields, one may use a 
single normalized complex three- vector, 



e(x) 



xi(x) 
x 2 {x) + iy 2 {x) 
^x 3 (x) +iy 3 (x) 



(8) 



where x\(x) > 0. For U(2), one may use two 
orthonormal complex five-vectors[^j. 

We have successfully used this method to sim- 
ulate U(l) and are now applying it to SU{2). 

5. FERMIONS WITHOUT DOUBLERS 

On the lattice each species of fermion appears 
as 16 different fcrmions. The root of this prob- 
lem is that the natural discretization of the Fermi 
action approximates the derivative by means of a 
gap of two lattice spacings. 

At the price of some nonlocality, we may leave 
out the unwanted states from the start B. Thus 
on a lattice of even size N = 2F, we may place 
independent fermionic variables , 0(2n) and ip{2n) 
only on the F 4 even sites 2n where n is a four- 
vector of integers, 



2n = (2ni, 2n 2 , 2n 3 , 2n^) 



(9) 
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and 1 < n, < F for i — 1, ... 4. To extend the vari- 
ables tp(2n) and ip(2n) to the nearest-neighbor 
sites 2n ± /t, we first define the Fourier variables 

4>{k) (and i)(k)), 



-y 



exp 



. 2im ■ k 



F 



ip{2n) 



(10) 



in which the sum extends over each rii , i = 
1, . . . , 4 from 1 to F and k is a four- vector of in- 
tegers k — (fci, k2, k^, k^) with 1 < ki < F. In 
terms of these Fourier variables, the dependent, 
nearest- neighbor variable ?p(2n + /}) is 



</-(2™ + A) = -LE' 



,iir(2n+fi)-k/F 



(11) 



in which the sum extends over each ki, i = 1, .., 4 
from 1 to F. In terms of the lattice delta functions 

F 



^(2j) = ^exp 



i=i 



F 



(12) 



we may write the nearest-neighbor variables as 
the sums 



F-r. 



ip{2n±fi) 



S ± (2j)^(2n + 2jfi). (13) 

j=l— n„ 



The action with independent fields at only F 4 
sites is 

5 = (2a) 4 ^^(271)1 -mip(2n) 



~ E ^ t^( 2n + A) - tf(2« - A)] },(14) 

in which the nearest- neighbor variables ip(2n±jl) 
are given by eq. ( [l3|) and the sum is over 1 < n,; < 
F as in eq.(|o|). 

We may now verify that there are no doublers. 
The Fourier series ([l0]) and ([ll]) diagonalize the 
action (114) and the lattice propagator is 



-ma + Y.^ln sin (^w) 



(15) 



For m = this propagator has poles only at k^ = 
F, which is the same point as fc M = 0. There are 



no doublers. Gauge fields may be added in the 
usual way or in the vectorial manner of section 
4. In a gauge theory with a gauge field U^n) 
defined on the link (ra, n + jl), one may construct 
the ordered product U(2n, 2n + 2j(i) of Wilson 
links U^n) from site 2n + 2jjl to site 2n for j > 
0. Thus by inserting the line U(2n,2n + 2jfi) 
into the action (p^), one may render it covariant. 
This procedure should also work for chiral gauge 
theories. 

Because of the lack of locality, the fermion 
propagator is not as sparse as the usual propa- 
gator. On the other hand, there are only one- 
sixteenth as many fermionic variables, and so the 
fermion propagator is smaller by a factor of 256. 
The present formalism of thinned fermions there- 
fore may be useful in practice as well as in prin- 
ciple. 

We intend to test this idea by simulating QED 
in two dimensions. 
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